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BRIDGE  HYDRAULIC  ANALYSIS  WITH  HEC-RAS‘ 


Vernon  Bonner  and  Gary  Brunner 
Hydrologic  Engineering  Center,  US  Army  Corps  of  Engineers 

Introduction 

The  Hydrologic  Engineering  Center  (HEC)  is  developing  next  generation  software  for 
one-dimensional  river  hydraulics.  The  HEC-RAS  River  Analysis  System  (EEC,  1995a)  is 
intended  to  be  the  successor  to  the  current  steady-flow  HEC-2  Water  Surface  Profiles  Program 
(HEC,  1990).  It  will  also  provide  unsteady  flow,  sediment  transport,  and  hydraulic  design 
capabilities  in  the  future.  A  single  data  definition  of  the  river  reach  and  bridge  data  is  used  for  all 
modeling  methods.  The  form  of  the  bridge  data  requirements,  the  hydraulic  routines,  and  the 
subsequent  research  into  bridge  flow  characteristics  and  modeling  approach  represent  a 
significant  improvement  in  one-dimensional  bridge  modeling.  This  paper  provides  an  overview 
of  HEC-RAS  representation  of  bridge  hydraulics  and  the  results  of  several  bridge  hydraulics 
investigations. 


HEC-RAS  Overview 

The  HEC-RAS  is  an  integrated  package,  designed  for  interactive  use  in  a  multi-tasking 
environment.  The  system  uses  a  Graphical  User  Interface  (GUI)  for  file  management,  data  entry 
and  editing,  program  execution,  and  output  display.  The  system  is  designed  to  provide  one¬ 
dimensional  river  modehng  using  steady-flow,  unsteady-flow  and  sediment-transport 
computations  based  on  a  single  geometric  representation  of  the  river  network.  Version  1 . 1 
provides  steady-flow  water  surface  profile  calculations  for  a  river  network  with  sub-critical, 
supercritical,  or  mixed-flow  regime 

The  program  has  been  developed  based  on  a  single  definition  of  the  river  geometric  data 
for  all  modeling.  The  five  steps  for  developing  a  hydraulic  model  are:  1)  create  a  project  file;  2) 
develop  the  river  network  and  enter  geometric  data;  3)  define  flow  and  boundary  conditions;  4) 
perform  hydraulic  analyses;  and  5)  review  results  and  produce  reports.  River  networks  are 
defined  by  drawing,  with  a  mouse,  a  schematic  of  the  river  reaches  from  upstream  to 
downstream.  As  reaches  are  connected  together,  junctions  are  automatically  formed  by  the 
program.  After  the  network  is  defined,  reach  and  junction  input  data  can  be  entered.  The  data 
editors  can  be  called  by  pressing  the  appropriate  icons  in  the  Geometric  Data  Window;  or  reach 
data  can  be  imported  from  HEC-2  data  sets. 

Cross-section  data  are  defined  by  reach  name  and  river  station.  Data  are  defined  by 
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station-elevation  coordinates,  up  to  500  coordinates  are  allowed.  There  is  no  maximum  number 
of  cross  sections.  The  cross  sections  are  stored  in  a  downstream  order  based  on  their  river-station 
number.  Cross  sections  can  be  easily  added  or  modified  in  any  order.  Cut,  copy,  and  paste 
features  are  provided,  along  with  separate  expansion  or  contraction  of  the  cross-section  elements 
of  overbanks  and  channel.  Cross-section  interpolation  is  provided  using  cross-section 
coordinates.  The  program  connects  adjacent  cross  sections  with  chords  at  the  boundaries,  bank 
stations,  and  minimum  elevation  point.  The  user  can  add  chords  graphically.  The  interpolated 
sections  are  marked  in  all  displays  to  differentiate  them  from  input  data. 

Steady-flow  data  are  defined  for  the  reach  at  any  cross-section  location.  Multiple-profile 
calculations  can  be  performed.  The  boundary  conditions  are  defined  at  downstream,  and/or 
upstream  ends  of  reaches  depending  on  flow  regime.  Internal  boundary  conditions  are  defined  at 
the  junctions.  Options  for  starting  profile  calculations  include:  known  water-surface  elevation, 
energy  slope  (normal  depth),  rating  curve,  and  critical  depth. 

Profile  calculations  are  performed  using  the  standard-step  procedure.  Overbank 
conveyance  is  computed  incrementally  at  coordinate  points  (HEC-2  style)  or  at  breaks  in 
roughness  (HEC-RAS  default).^  Subcritical,  supercritical,  and  mixed-flow  profile  calculations 
can  be  performed.  The  critical-depth  routine  searches  the  entire  range  of  depths  and  locates 
multiple  specific  energy  minima.  The  location  of  the  transition  between  supercritical  and 
subcritical  flow  is  determined  based  on  momentum  calculations.  Detailed  hydraulic  jump 
location  and  losses  are  not  computed;  however,  the  jump  location  is  defined  between  two 
adjacent  cross  sections. 

Tabular  output  is  available  using  pre-defined  and  user-defined  tables.  Cross-section 
tables  provide  detailed  hydraulic  information  at  a  single  location,  for  a  profile.  Profile  tables 
provide  summary  information  for  all  locations  and  profiles.  Pre-defined  summary  tables  are 
available  for  the  cross  section,  bridge,  and  culvert  computations.  User-defined  tables  can  be 
developed,  from  a  menu  of  170  output  variables,  and  stored  for  future  use  like  the  pre-defined 
tables. 


Graphical  displays  are  available  for  cross  sections,  profiles,  rating  curves,  and  a  X-Y-Z 
perspective  plot  of  the  river  reach,  as  shown  in  Figure  1.  The  geometric  data  and  computed 
results  can  be  displayed  from  the  View  option,  provided  in  most  of  the  data-input  editors.  User 
control  is  provided  for  variables  to  plot,  line  color,  width  and  type,  plus  axis  labels  and  scales. 
The  user  can  also  zoom-in  on  selected  portions  of  the  display,  and  zoom-out  to  the  original  size. 
All  graphics  are  in  vector  form  using  calls  to  the  Window's™  Graphics  Device  Interface. 
Graphics  can  be  sent  to  output  devices  through  the  Windows  print  manager,  or  they  can  be 
written  to  a  meta  file  or  sent  to  the  Windows  chp  board. 


^  Comparisons  with  HEC-2  results  and  HEC-RAS  using  both  conveyance  calculation 
methods  were  presented  atASCE  Hydraulic  Engineering  ‘94  (Bonner  &  Brunner,  1994). 
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Bridge  and  Culvert  Routines 

The  bridge  routines  in  HEC-RAS  enables  analysis  of  bridge  hydraulics  by  several 
different  methods  without  changing  the  bridge  geometry.  The  model  utilizes  four  user-  defined 
cross  sections  in  the  computations  of  energy  losses  due  to  the  stmcture,  as  shown  in  Figure  2. 

An  effective-area  option  is  used  with  the  bounding  cross  sections  2  and  3  to  define  the  ineffective 
flow  areas  in  those  sections  (zones  marked  A-B  and  C-D  in  Figure  2).  Cross  sections  are 
formulated  inside  the  bridge  by  combining  bridge  geometry  with  the  two  bounding  cross 
sections.  The  bridge  geometry  is  defined  by  the  roadway/deck,  piers,  and  abutments  as  separate 
input  items.  Bridge  data  are  entered  through  the  editor,  shown  in  Figure  3.  The  bridge  modeling 
methods  are  also  selected  from  a  memu  of  options,  called  from  the  bridge  editor. 

Low-Flow  Computations.  The  program  first  uses  the  momentum  equation  to  define  the 
class  of  flow.  For  Class  A  low  flow  (completely  subcritical),  the  modeler  can  select  any  or  all  of 
the  following  three  methods  to  compute  bridge  energy  losses:  standard-step  energy;  momentum; 
or  Yamell  equation.  The  USGS-Federal  Highway  WSPRO  bridge  routine  (FHWA,  1990)  will  be 
included  in  a  later  program  release.  If  more  than  one  method  is  selected,  the  user  must  choose  a 
single  method  as  the  final  solution,  or  select  the  method  that  produces  the  highest  energy  loss 
through  the  structure.  For  Class  B  low  flow  (passes  through  critical  depth)  the  program  uses  the 
momentum  equation.  Class  C  low  flow  (completely  supercritical)  can  be  modeled  with  either  the 
standard  step  energy  method  or  the  momentum  equation.  Also,  the  program  can  incorporate 
forces  due  to  water  weight  and/or  friction  components  in  addition  to  the  pier  impact  losses  for 
momentum. 
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Figure  2.  Cross  Sections  to  Model  Bridge  Flow  Transitions 

Pressure  Flow.  When  the  flow  comes  into  contact  with  the  low  cord  of  the  bridge, 
pressure  flow  begins.  The  program  uses  energy-based  (like  HEC-2  Normal  Bridge)  or  pressure- 
flow  equations.  It  checks  for  the  possibility  of  pressure  flow  when  the  upstream  energy-grade 
line  goes  above  the  maximum  low  chord.  The  user  has  the  option  of  using  the  water  surface 
instead  of  energy  elevation.  The  program  will  handle  two  cases  of  pressure  flow.  When  only  the 
upstream  side  of  the  bridge  is  in  contact  with  the  flow,  the  sluice  gate  equation  is  used.  If  both 
the  upstream  and  downstream  sides  of  the  bridge  are  in  contact  with  the  flow,  the  full-flow 
orifice  equation  is  used. 

Weir  Flow.  When  water  flows  over  the  bridge  and/or  roadway,  the  overflow  is  calculated 
using  a  standard  weir  equation.  For  high  tailwater  conditions,  the  amount  of  weir  flow  is  reduced 
to  account  for  the  effects  of  submergence.  If  the  weir  becomes  highly  submerged,  the  program 
will  switch  to  calculating  energy  losses  by  the  standard-step  energy  method.  The  criterion  for 
switching  to  energy-based  calculations  is  user  controllable.  When  combinations  of  low  flow  or 
pressure  flow  occur  with  weir  flow,  an  iterative  procedure  is  used  to  determine  the  amount  of 
each  type  of  flow. 
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Culvert  hydraulics.  The 
modeling  approach  for  culverts  is 
similar  to  that  for  bridges.  The 
cross-section  layout,  the  use  of 
ineffective  areas,  and  the  selection 
of  contraction  and  expansion 
coefficients  are  the  same.  For 
culvert  hydraulics,  the  program 
uses  the  Federal  Highway 
Administrations  (FHWA,  1985) 
culvert  equations  to  model  inlet 
control.  Outlet  control  is  analyzed 
by  either  direct-step  backwater 
calculations  or  full-flow  friction 
losses,  plus  entrance  and  exit 
losses.  The  culvert  routines  have 
the  ability  to  model  the  following 
shapes:  box;  circular;  arch;  pipe 
arch;  and  elliptical.  Multiple 
culverts,  of  different  types,  can  be 
modeled  for  a  single  location. 


Figure  3.  HEC-RAS  Bridge  Data  Editor 


Multiple  bridge  openings.  Multiple  openings  can  be  modeled  by  two  approaches,  as 
divided  flow  in  two  reaches  or  by  the  multiple-opening  approach.  The  multiple-opening 
approach  can  analyze  combinations  of  three  types  of  openings:  bridges,  culvert  groups,  and 
conveyance  areas.  Up  to  seven  openings  can  be  defined  at  any  one  river  crossing.  Each  opening 
is  evaluated  separately  and  the  total  flow  is  distributed  such  that  the  energy  loss  in  each  is  equal. 


Bridge  Flow  Transitions 


A  Master  of  Science  thesis  project  (Hunt,  1995)  was  conducted  to  investigate  bridge 
expansion  and  contraction  reach  lengths  and  coefficients.  Two-dimensional  models  of  ideahzed 
bridge  crossings  were  developed  using  RMA-2V  (King,  1994).  River  slopes,  bridge  opening 
widths,  overbank  to  channel  n-value  ratios,  and  abutment  type  were  varied;  a  total  of  76  cases 
were  modeled.  From  the  model  results,  regression  analyses  were  performed  to  develop  predictor 
equations  for  contraction  and  expansion  reach  lengths,  ratios,  and  coefficients.  (Figure  4 
illustrates  the  transition  lengths  and  ratios.)  The  results  and  general  recommendations  from  this 
modeling  review  are  summarized  here  from  HEC  Research  Document  No.  42  (HEC,  1995). 
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typical  flow 
transition  pattern 


Flow 


assumed  flow  transition 
pattern  for  1 -dimensional 


Figure  4.  Conceptual  Illustration  of  Transition  Reaches 


The  flow  transitions  through  a  bridge  crossing,  that  blocks  a  portion  of  the  overbank  area, 
is  typically  modeled  with  four  cross  sections.  The  downstream  and  upstream  sections  1  and  4 
represent  the  full  floodplain  conveyance.  The  bridge-bounding  cross  sections  2  and  3  represent 
the  effective  flow  area  just  downstream  and  upstream  from  the  bridge.  The  bridge  is  modeled 
with  the  bounding  cross  sections  and  additional  bridge  data.  The  question  is  how  far  does  it  take 
for  the  flow  transition  to  occur,  and  where  to  locate  the  full-flow  cross  sections  1  and  4. 

Expansion  Reach  lengths  (LJ.  The  expansion  ratio  (ER  in  Fig.  4)  was  less  than  4: 1  for 
all  of  the  idealized  cases  modeled.  The  mean  and  median  values  of  the  expansion  ratio  for  the 
idealized  cases  were  both  around  1.5:1.  The  idealized  cases  included  a  wide  range  of  hydraulic 
and  geometric  conditions.  These  observations  indicate  that  the  traditional  4: 1  rule  of  thumb  will 
over  predict  the  expansion  reach  length  for  most  situations. 

Many  independent  variables  and  combinations  of  variables  were  investigated  to  find  a 
possible  correlation  with  L^.  The  variable  which  showed  the  greatest  correlation  was  the  ratio  of 
the  main  channel  Froude  number  at  the  most  constricted  Section  2  to  that  at  the  normal  flow 
Section  1.  (See  Figure  4  for  cross-section  references.)  The  best-fitting  equation  for  the 
expansion  reach  length  is: 
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(1) 


Lg  =  -  298  +  257  (-^)  +  0.918  (L^^s)  +  0.00479  (Q) 

^Cl 

for  which  =0.84  and  =  96  feet,  with 

Lg  =  length  of  the  expansion  reach,  in  feet, 

Fc2  =  main  channel  Froude  number  at  Section  2, 

Fci  =  main  channel  Froude  number  at  Section  1, 

=  average  length  of  bridge  obstruction,  in  feet, 

Q  =  total  discharge,  cfs, 

=  the  adjusted  determination  coefficient  (the  percentage  of  variance  of  the 
dependent  variable  from  the  mean  which  is  explained  by  the  regression 
equation),  and 

Sg  =  standard  error  of  estimate. 


Similarly,  the  regression  equation  for  the  expansion  ratio  was  found  to  be 


ER 


^  =  0.421  +  0.485  (-^)  +  1.80x10  ®  (Q) 

^obs 


(2) 


for  which  R^  =  0.71  and  Sg  =  0.26. 

Figures  5  and  6  are  plots  of  the  observed  values  versus  those  predicted  by  the  regression 
equations  for  and  ER,  respectively. 
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In  some  studies,  a  high  level  of  sophistication  in  the  evaluation  of  the  transition  reach 
lengths  is  not  justified.  For  such  studies,  and  for  a  starting  point  in  more  detailed  studies.  Table 
1  offers  ranges  of  expansion  ratios  which  can  be  used  for  different  degrees  of  constriction, 
different  slopes,  and  different  ratios  of  overbank  roughness  to  main  channel  roughness.  Once 
an  expansion  ratio  is  selected,  the  distance  to  the  downstream  end  of  the  expansion  reach  (the 
location  of  Section  1  on  Figure  4)  is  found  by  multiplying  the  expansion  ratio  by  the  average 
obstruction  length.  The  average  obstruction  length  is  half  of  the  total  reduction  in  floodplain 
width  caused  by  the  two  bridge  approach  embankments. 

In  Table  1,  b/B  is  the  ratio  of  the  bridge  opening  width  to  the  total  floodplain  width, 
is  the  Manning  n  value  for  the  overbank,  is  the  n  value  for  the  main  channel,  and  S  is  the 
longitudinal  slope.  The  values  in  the  interior  of  the  table  are  the  ranges  of  the  expansion  ratio. 
For  each  range,  the  higher  value  is  typically  associated  with  a  higher  discharge. 

Table  1.  Ranges  of  Expansion  Ratios 


b/B  =  0.10 

Slope: 

«ob  /  «c  =  1 

Mob  /  =  2 

Mob/ Me  =  4 

1  fi^mile 

1.4 -3.6 

1.3 -3.0 

1.2 -2.1 

5  ft/mile 

1.0 -2.5 

0.8  -  2.0 

0.8  -  2.0 

1 0  ft/mile 

1.0 -2.2 

0.8 -2.0 

0.8  -  2.0 

b/B  =  0.25 

1  ft/mile 

1.6 -3.0 

1.4 -2.5 

1.2 -2.0 

5  ft/mile 

1.5 -2.5 

1.3 -2.0 

1.3 -2.0 

10  ft/mile 

1.5 -2.0 

1.3  -2.0 

1.3 -2.0 

b/B  =  0.50 

1  fl/mile 

1.4 -2.6 

1.3  -  1.9 

1.2- 1.4 

5  fl/mile 

1.3 -2.1 

1.2- 1.6 

1.0- 1.4 

10  ft/ mile 

1.3 -2.0 

1.2- 1.5 

1.0 -1.4 

The  ranges  in  Table  1 ,  as  well  as  the  ranges  of  other  parameters  to  be  presented  later,  capture  the 
ranges  of  the  idealized  model  data  from  this  study.  Extrapolation  of  expansion  ratios  for 
constriction  ratios,  slopes  or  roughness  ratios  outside  of  the  ranges  used  in  this  table  should  be 
done  with  care.  The  expansion  ratio  should  not  exceed  4:1,  nor  should  it  be  less  than  0.5:1, 
unless  there  is  site-specific  field  information  to  substantiate  such  values.  The  ratio  of  overbank 
roughness  to  main-channel  roughness  provides  information  about  the  relative  conveyances  of  the 
overbank  and  main  channel. 
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The  user  should  note  that  in  the  data  used  to  develop  these  recommendations,  all  cases 
had  a  main-channel  n  value  of  0.04.  For  significantly  higher  or  lower  main-channel  n  values,  the 
n  value  ratios  will  have  a  different  meaning  with  respect  to  overbank  roughness.  It  is  impossible 
to  determine  from  the  data  of  this  study  whether  this  would  introduce  significant  error  in  the  use 
of  these  recommendations. 

It  is  recommended  that  the  user  start  with  an  expansion  ratio  from  Table  1,  locate  Section 
1  according  to  that  expansion  ratio,  set  the  main  channel  and  overbank  reach  lengths  as 
appropriate,  and  limit  the  effective  flow  area  at  Section  2  to  the  approximate  bridge  opening 
width.  The  program  should  then  be  run  and  the  main  channel  Froude  numbers  at  Sections  2  and 
1  read  from  the  model  output.  Use  these  Froude  number  values  to  determine  a  new  expansion 
length  from  the  appropriate  equation,  move  Section  1  as  appropriate  and  recompute.  Unless  the 
geometry  is  changing  rapidly  in  the  vicinity  of  Section  1,  no  more  than  two  iterations  after  the 
initial  run  should  be  required. 

When  the  expansion  ratio  is  large,  say  greater  than  3:1,  the  resulting  reach  length  may  be 
so  long  as  to  require  intermediate  cross  sections  which  reflect  the  changing  width  of  the  effective 
flow  area.  These  intermediate  sections  are  necessary  to  reduce  the  reach  lengths  when  they 
would  otherwise  be  top  long  for  the  linear  approximation  of  energy  loss  that  is  incorporated  in 
the  standard  step  method.  These  interpolated  sections  are  easy  to  create  in  the  HEC-RAS 
program,  because  it  has  a  graphical  cross  section  interpolation  feature.  The  importance  of 
interpolated  sections  in  a  given  reach  can  be  tested  by  first  inserting  one  interpolated  section  and 
seeing  the  effect  on  the  results.  If  the  effect  is  significant,  the  subreaches  should  be  subdivided 
into  smaller  units  until  the  effect  of  further  subdivision  is  inconsequential. 

Contraction  Reach  Lengths  (L^.).  In  contrast  to  the  expansion  reach  length  results,  the 
results  for  contraction  lend  some  support  to  the  traditional  rule  of  thumb  which  recommends  the 
use  of  a  1:1  contraction  ratio.  The  range  of  values  for  this  ratio  was  from  0.7: 1  to  2.3:1.  The 
median  and  mean  values  were  both  around  1.1  to  1. 

The  Froude  number  ratio  in  the  previous  two  equations  also  proved  to  be  significant  in  its 
relationship  to  the  contraction  reach  length.  Surprisingly,  the  Froude  number  ratio  which 
involved  the  upstream  (Section  4)  Froude  number  did  not  have  as  strong  a  correlation  as  the  one 
involving  the  Section  1  value.  Here  again  the  degree  of  constriction,  in  comparison  with  the 
undisturbed  flow  condition,  is  of  high  significance.  The  most  significant  independent  variable 
for  this  parameter,  however,  was  the  percentage  of  the  total  discharge  conveyed  by  the  two 
overbanks.  The  best-fit  equation  from  the  regression  analysis  is 

L  =  263 +38.8(^^) +257(^)2-58.7  (-^f  5 +0.161  (3) 

with  R^  =  0.87  and  S^.  =  3 1  feet.  In  this  equation 
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Qob  =  the  discharge  conveyed  by  the  two  overbank  sections,  in  cfs,  and 
=  the  Manning  n  value  for  the  overbank  sections. 

Figure  7  shows  the  observed  versus  predicted  values  for  Equation  3.  The  contraction 
length  values  did  not  vary  much  as  a  function  of  the  bridge  opening  width  or  the  average 
obstruction  length.  As  a  result  most  of  the  cases  with  the  widest  opening  width,  and  therefore  the 
shortest  average  obstruction  length,  had  the  highest  contraction  ratios.  The  numerator  of  the  ratio 
varied  only  slightly  while  the  denominator  varied  greatly. 


Figure  7.  Goodness-of-Fit  Plot  for  Contraction  Reach  Length  Regression 
Equation  (Equation  3). 


None  of  the  attempted  regression  relationships  were  good  predictors  of  the  contraction 
ratio.  Equation  4  provided  the  best  fit  of  all  the  combinations  of  independent  variables  tried. 
The  scatter  around  a  CR  value  of  1  is  shown  in  Figure  8,  a  plot  of  the  observed  versus  predicted 
values  of  the  contraction  ratio.  The  regression  equation  for  the  contraction  ratio  is: 

F"  Q  n 

CR  =  1.4  -  0.333(-^)  +  1.86(-^)2  -  0.19(-5^)°-s  (4) 

F,/  'q' 


This  equation  has  an  =  0.65  and  =0.19. 

An  unfortunate  feature  of  this  equation  is  the  negative  sign  on  the  Froude  number  ratio 
term.  This  negative  term  indicates  that  the  contraction  ratio  should  become  smaller  as  the 
constriction  gets  more  severe.  Given  the  low  R^  for  equation  4,  equation  3  is  preferable. 
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Figure  8.  Goodness-of-Fit  Plot  for  Contraction  Ratio  Regression  Equation 
(Equation  4). 


Ranges  of  contraction  ratios  for  different  conditions  are  presented  in  Table  2.  An  average 
value  or  a  value  from  the  table  can  be  used  as  a  starting  estimate  and  contraction  reach  length  f  or 
studies  which  do  not  justify  a  more  detailed  evaluation. 


Table  2.  Ranges  of  Contraction  Ratios 


Wob  /  =  1 

IBHBiSIIIH 

”ob/«c=  4 

S  =  1  ft/mile 

1.0  -  2.3 

00 

d 

0.7  -  1.3 

5  fl/mile 

1.0- 1.9 

0.8  -  1.5 

0.7  -  1.2 

10  ft/mile 

1.0-  1.9 

0.8  -  1.4 

0.7  -  1.2 

When  the  conditions  are  within  or  near  those  of  the  data,  the  contraction  reach  length 
regression  Equation  3,  may  be  used.  In  cases  where  the  floodplain  scale  and  discharge  are 
significantly  larger  or  smaller  than  those  that  were  used  in  developing  the  regression  formulae, 
Equation  3  should  not  be  used.  The  recommended  approach  for  estimating  the  contraction  ratio 
at  this  time  is  to  compute  a  value  from  Equation  4  and  check  it  against  the  values  in  Table  2.  As 
with  the  expansion  reach  lengths,  the  modeler  must  use  these  equations  and  the  values  from 
Table  2  with  extreme  caution  when  the  prototype  is  outside  of  the  range  of  data  used  in  this 
study.  The  contraction  ratio  should  not  exceed  2.5:1  nor  should  it  be  less  than  0.3:1. 

FYpansinn  Coefficients  (CJ.  Unlike  the  transition  reach  lengths,  the  transition 
coefficients  did  not  lend  themselves  to  strong  regression  relationships.  This  situation  is  partly 
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due  to  the  fact  that  the  velocity  head  differences  were  so  small  in  many  instances  as  to  render  the 
coefficient  values  insignificant.  Calibration  of  the  coefficients  under  these  conditions  is 
obviously  meaningless.  Despite  these  difficulties,  some  trends  were  apparent  in  the  expansion 
coefficient.  The  ratio  of  the  hydraulic  depth  on  the  overbanks  to  the  hydraulic  depth  in  the  main 
channel  showed  some  correlation  with  Cg.  The  best  regression  relationship  was: 

Cg  =  -0.092  +  0.570  (-?^)  +  0.075  (-^)  (5) 

^Cl 


for  which  =  0.55  and  Se  =  0.10,  with 

Dob  =  hydrauUc  depth  (flow  area  divided  by  top  width)  for  the  overbank  at  the 

normal  flow  section  (Section  1),  in  feet,  and 

Dg  =  hydraulic  depth  for  the  main  channel  at  the  normal  flow  section  (Section 

1),  in  feet. 

Figure  9  shows  the  goodness  of  fit  for  this  equation.  The  analysis  of  the  data  with  regard  to  the 
expansion  coefficients  did  not  yield  a  regression  equation  which  fit  the  data  well.  Equation  5 
was  the  best  equation  obtained  for  predicting  the  value  of  this  coefficient.  The  calibrated 
expansion  coefficients  ranged  from  0.1  to  0.65.  The  median  value  was  0.3.  Recalling  that  the 
traditional  rule  of  thumb  for  this  coefficient  suggests  a  standard  value  of  0.5,  it  appears  that  the 
application  of  this  rule  could  lead  to  an  over  prediction  of  energy  loss  in  the  expansion  reach. 


Figure  9.  Goodness-of-Fit  Plot  for  Expansion  Coefficient  Regression  Equation 
(Equation  5). 
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It  is  recommended  that  the  modeler  use  an  average  value  or  Equation  5  to  find  an  initial 
value,  then  perform  a  sensitivity  analysis  using  values  of  the  coefficient  that  are  0.2  higher  and 
0.2  lower  than  the  initial  value.  The  plus  or  minus  0.2  range  defines  the  95%  confidence  band 
for  Equation  5  as  a  predictor  within  the  domain  of  the  regression  data.  If  the  difference  in  results 
between  the  two  ends  of  this  range  is  substantial,  then  the  conservative  value  should  be  used. 

The  expansion  coefficient  should  higher  than  the  contraction  coefficient  and  less  than  1.0. 

Contraction  Coefficients  (C^) .  Of  the  76  cases  used  in  the  regression  analysis  (those  with 
symmetric  openings  and  spUl-through  abutments),  69  had  calibrated  C^  values  of  0.10.  These 
included  cases  for  which  the  contraction  coefficient  had  no  appreciable  significance,  as  well  as 
the  28  cases  wherein  the  RMA-2  water  surface  elevation  at  the  approach  section  was  too  low  to 
be  reached  in  HEC-RAS.  Because  of  these  conditions,  the  regression  analysis  was  unfruitful.  In 
addition  to  the  regression  study  with  all  of  the  data,  an  attempt  at  regression  was  made  which 
incorporated  only  20  cases.  For  this  analysis  those  cases  in  which  the  contraction  coefficient  was 
inconsequential  were  omitted.  This  attempt  also  failed  to  yield  a  satisfactory  regression 
relationship  because  13  of  these  20  cases  still  had  calibrated  coefficient  values  of  0.10. 

The  values  for  the  contraction  coefficient  ranged  from  0. 10  to  0.50.  The  mean  was  0. 12 
and  the  median  value  obviously  was  0. 10.  Here  again  is  a  suggestion  that  the  traditional 
standard  value  used  for  bridges,  0.30,  is  probably  too  high. 

The  data  of  this  study  did  not  lend  itself  to  regression  of  the  contraction  coefficient 
values.  For  nearly  all  of  the  cases  the  value  that  was  determined  was  0.1,  which  was  considered 
to  be  the  minimiiTn  acceptable  value.  The  following  table  presents  recommended  ranges  of  the 
contraction  coefficient  for  various  degrees  of  constriction,  for  use  in  the  absence  of  calibration 
information. 


Table  3.  Contraction  Coefficient  Values 


Degree  of  Constriction 

Recommended  Contraction  Coefficient 

0%  <  b/B  <  25% 

0.3  -  0.5 

25%<b/B<50% 

0.1 -0.3 

50%  <  b/B  <100% 

0.1 

■Sjiimmarv  The  preceding  recommendations  represent  a  substantial  improvement  over  the 
guidance  information  that  was  previously  available  on  the  evaluation  of  transition  reach  lengths 
and  coefficients.  They  are  based  on  model  data  which,  like  all  data,  have  a  limited  scope  of 
direct  application.  Certain  situations,  such  as  highly  skewed  bridge  crossings  and  bridges  at 
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locations  of  sharp  curvature  in  the  floodplain  were  not  addressed  by  this  study.  Even  so,  these 
recommendations  may  be  applicable  to  such  situations  if  proper  care  is  taken  and  good 
engineering  judgement  is  employed. 

In  applying  these  recommendations,  the  modeler  should  always  consider  the  range  of 
hydraulic  and  geometric  conditions  included  in  the  data.  Wherever  possible,  the  transition  reach 
lengths  used  in  the  model  should  be  validated  by  field  observations  of  the  site  in  question, 
preferably  under  conditions  of  high  discharge.  The  evaluation  of  contraction  and  expansion 
coefficients  should  ideally  be  substantiated  by  site-specific  calibration  data,  such  as  stage- 
discharge  measurements  upstream  of  the  bridge.  The  guidelines  given  here  recognize  the  fact 
that  site-specific  field  information  is  often  unavailable  or  very  expensive  to  obtain. 

Program  Testing 

Initial  program  testing  compared  HEC-RAS  with  the  results  from  the  current  HEC-2 
program,  primarily  evaluating  differences  in  cross-section  water-surface  elevations.  These  tests 
only  involved  cross-sectional  data;  no  bridges  or  culverts  were  in  the  models.  Water  surface 
profile  comparisons  were  made  for  both  the  HEC-2  and  HEC-RAS  approaches  for  conveyance 
calculations  (Brunner  &  Bonner,  1994).  Ninety-six  percent  of  HEC-RAS  elevations  were  within 
6  mm  of  HEC-2  profile  results,  when  using  the  HEC-2  conveyance  computation  option. 

However,  the  computed  profiles  are  higher,  than  those  from  HEC-2,  when  the  default 
conveyance  computation  is  used.  We  believe  that  the  default  HEC-RAS  method  is  more 
consistent  with  the  theory  and  with  computer  programs  HEC-6  (HEC,  1993),  UNET  (HEC, 
1995d),  and  WSPRO  (FHWA,  1990). 

The  bridge  routines  of  HEC-RAS,  HEC-2,  and  WSPRO  were  tested  using  21  USGS  data 
sets  from  the  Bay  St.  Louis  Laboratory  (Ming  et  al.,  1978).  “In  general,  all  models  were  able  to 
compute  water  surface  profiles  within  the  tolerance  of  the  observed  data”  (HEC,  1995c).  “The 
variation  of  the  water  surface  at  any  given  cross  section  was  on  the  order  of  0.1  to  0.3  feet” 

(ibid).  For  HEC-RAS  and  HEC-2,  the  energy-based  methods  reproduced  observed  bridge  low- 
flow  losses  better  than  the  Yamell  method.  Also,  the  apparent  downstream  expansion  reach- 
lengths  were  shorter  than  rates  suggested  in  HEC  guidehnes  (HEC,  1990).  Because  all  the 
prototype  bridge  data  came  from  similar  wide,  heavily-vegetated  flood  plains  with  low  velocities, 
additional  research  was  conducted  using  the  RMA-2V  computer  program  (King,  1994). 
Application  of  the  2-D  model  to  the  prototype  data  demonstrated  that  RMA-2V  could  reproduce 
observed  bridge-flow  depths  and  transitions.  Based  on  these  results,  the  2-D  program  was  used 
to  gimnlatp.  different  bridge  configurations,  as  reported  in  the  section  Bridge  Flow  Transitions. 

Along  with  internal  HEC  testing,  two  Beta  versions  were  released  to  the  public  for  testing 
(250  for  Beta  1  and  200  for  Beta  2).  After  all  the  testing  was  completed,  and  final  corrections 
and  additions  were  made,  HEC  released  Version  1.0  in  July  1995.  Program  development 
continued,  as  well  as  error  corrections  found  in  Version  1.0.  The  program  update.  Version  1.1, 
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was  released  in  January  1996.  Program  development  on  the  steady-flow  model  is  expected 
through  fiscal  year  1996,  when  most  of  the  expected  features  will  be  completed.  After  that,  work 
will  proceed  with  the  unsteady-flow  model. 
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